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Abstract. We present and study lattice and off- lattice microscopic models in which 
particles interact via a local anisotropic rule. The rule induces preferential hopping 
along one direction, so that a net current sets in if allowed by boundary conditions. 
This may be viewed as an oversimplification of the situation concerning certain traffic 
and fiow problems. The emphasis in our study is on the infiuence of dynamic details on 
the resulting (non-equilibrium) steady state. In particular, we shall discuss on the simi- 
larities and differences between a lattice model and its continuous counterpart, namely, 
a Lennard-Jones analogue in which the particles' coordinates vary continuously. Our 
study, which involves a large series of computer simulations, in particular reveals that 
spatial discretization will often modify the resulting morphological properties and even 
induce a different phase diagram and criticality. 

1 Introduction 

Many systems out of equilibrium |1I2| exhibit spatial striped patterns on macro- 
scopic scales. These are often caused by transport of matter or charge induced 
by a drive which leads to heterogeneous ordering. Such phenomenology occurs in 
flowing fluids 3 , and during phase separation in colloidal II , granular |5I6| , and 
liquid-liquid mixtures. Further examples are wind ripples in sand trails by 
animals and pedestrians and the anisotropics observed in high temperature 
superconductors 'lU'll' and in two-dimensional electron gases |12I13| . 

Studies of these situations, often described as nonequilibrium phase tran- 
sitions, have generally focused on lattice systems |14I15I16I17I18| . i.e., models 
based on a discretization of space and in considering interacting particles that 
move according to simple local rules. Such simplicity sometimes allows for ex- 
act calculations and is easy to be implemented in a computer. Moreover, some 
powerful techniques have been developed to deal with these situations, includ- 
ing nonequilibrium statistical field theory. However, lattice models are perhaps 
a too crude oversimplification of fluid systems so that the robustness of such an 
approach merits a detailed study. 

The present paper describes Monte Carlo (MC) simulations and field theoret- 
ical calculations that aim at illustrating how slight modifications of dynamics at 
the microscopic level may influence, even quantitatively, the resulting (nonequi- 
librium) steady state. We are also, in particular concerned with the influence of 
dynamics on criticality. With this objective, we take as a reference the driven 



2 M. Di'ez-Minguito et al. 



lattice gas (DLG), namely, a kinetic noncquilibrium Ising model with conserved 
dynamics. This system has become a prototype for anisotropic behavior, and 
it has been useful to model, for instance, ionic currents ^\ and traffic flows 
\W\ . In fact, in certain aspects, this model is more realistic for traffic flows than 
the standard asymmetric simple exclusion process TTO. Here we compare the 
transport and critical properties of the DLG with those for apparently close 
lattice and off-lattice models. There is some related previous work addressing 
the issue of how minor variations in the dynamics may induce dramatic mor- 
phological changes both in the early time kinetics and in the stationary state 
|2()I21I22] . However, these papers do not focus on transport nor on critical prop- 
erties. We here in particular investigate the question of how the lattice itself 
may condition transport, structural and critical properties and, with this aim, 
we consider nearest-neighbor (NN) and next-nearest-neighbor (NNN) interac- 
tions. We also compare with a microscopically off-lattice representation of the 
driven lattice gas in which the particles' spatial coordinates vary continuously. 
A principal conclusion is that spatial discretization may change significantly not 
only morphological and early-time kinetics properties, but also critical proper- 
ties. This is in contrast with the concept of universality in equilibrium systems, 
where critical properties are independent of dynamic details. 

2 Driven Lattice Gases 

The driven lattice gas, initially proposed by Katz, Lebowitz, and Spohn 23 , is a 
nonequilibrium extension of the Ising model with conserved dynamics. The DLG 
consists of a rf-dimensional square lattice gas in which pair of particles interact 
via an attractive and short-range Ising-like Hamiltonian, 

U-k) 

Here ak = 0(1) is the lattice occupation number at site k for an empty (occupied) 
state and the sum runs over all the NN sites (the accessible sites are depicted 
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Fig. 1. Schematic comparison of the sites a particle (at the center, marked with a dot) 
may occupy (if the corresponding site is empty) for nearest-neighbor (NN) and next- 
nearest-neighbor (NNN) hops at equilibrium (left) and in the presence of an "infinite" 
horizontal field (right). The particle-hole exchange between neighbors is either forbid- 
den (x) or allowed (^), depending on the field value. 
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in Fig.^. Dynamics is induced by the competion between a heat bath at tem- 
perature T and an external driving field E which favors particle hops along one 
of the principal lattice directions, say horizontally {x), as if the particles were 
positively charged. Consequently, for periodic boundary conditions, a nontrivial 
nonequilibrium steady state is set in asymptotically. MC simulations by a biased 
Metropolis rate reveal that, as in equilibrium, the DLG undergoes a second order 
phase transition. At high enough temperature, the system is in a disordered state 
while, below a critical point (at T < Te) it orders displaying anisotropic phase 
segregation. That is, an anisotropic (striped for d = 2) rich-particle phase then 
coexists with its gas. It is also found that the critical temperature Te monotoni- 
cally increases with E. More specifically, for d = 2, assuming a half filled square 
lattice in the large field limit (in order to maximize the nonequilibrium effect), 
one has a nonequilibrium critical point at Too ~ \ATq, where the equilibrium 
value is Tq — 2.269Jfc^^. It was numerically shown that this belongs to a univer- 
sality class other than the Onsager one, e.g., MC data indicates that the order 
parameter critical exponent is /3dlg — 1/3 |17I24| (instead of the Onsager value 
1/8). 

Other key features concern the two-particle correlation function C{x,y) and 
its Fourier transform S{kx, ky), i.e., the structure factor. As depicted in the left 
graph of Fig. |21 correlations are favored (inhibited) along (against) the field 
direction. In fact, the DLG shows a slow decay of the two-point correlations 
due to the spatial anisotropy associated with the dynamics [33] . This long range 
behavior translates into a characteristic discontinuity singularity at the origin 
(limfc^^o 'S'li 7^ limfe^^o '5'_l) in the structure factor ^16,. which is confirmed in 
Fig.EI 

How do all these features depend on the number of neighbor sites to which a 
particle can hop? Or in other words, how robust is the behavior when extending 
interactions and accessible sites to the NNN? 
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Fig. 2. Parallel (squares) and transverse (triangles) components of the two-point cor- 
relation function (left) and the structure factor (right) above criticality with NN (filled 
symbols) and NNN (empty symbols) interactions for a 128 x 128 half filled lattice. The 
inset shows the x^^ power law decay in for both discrete cases: DLG (o) and NDLG 
(X). 



4 M. Di'ez-Minguito et al. 



Previous work has shown that extending hopping in the DLG to NNN leads 
to an inversion of triangular anisotropics during the formation of clusters |21| . 
and also that dramatic changes occur in the steady state, including the fact that, 
contrary to the DLG with NN interactions, the critical temperature decreases 
with increasing E |22. However, other important features such as correlations 
and criticality seem to remain invariant. Analysis of the parallel {C'\\) and trans- 
verse (C_l) components reveals that correlations are quantitatively similar for the 
DLG and for the DLG with NNN interactions (henceforth NDLG) —although 
somehow weaker for the latter case. Also persists a slow decay of correlations 
which yield to the discontinuity at the origin of S{kx, ky). These facts are shown 
in Fig. 121 

On the other hand, recent MC simulations of the NDLG indicate that the 
order parameter critical exponent is /3ndlg ~ 1/3 26 , as for the DLG. The 
anisotropic diffusive system approach '28' , which is a Langevin-type (meso- 
scopic) description, predicts this critical behavior. In both cases, DLG and 
NDLG, the Langevin equations, as derived by coarse graining the master equa- 
tion, lead to (3 = 1/3. These two Langevin equations are identical, except for 
new entropic terms in the NDLG due to the presence of additional neighbors 

m 

The fact that extending particle hops and interaction to the diagonal sites 
leaves invariant both correlations and criticality seems to indicate that the two 
systems, DLG and NDLG, belong to the same universality class. 



3 A Driven Off-lattice Gas 

In order to deep further on this interesting issue, we studied to what extent 
the DLG behavior depends on the lattice itself. With this aim, we considered a 
driven system with continuous variation of the particles' spatial coordinates — 
instead of the discrete variations in the DLG — which follows as close as possible 
the DLG strategy. In particular, we analyzed an off-lattice, microscopically- 
continuum analog of the DLG with the symmetries and short-range interaction 
of this model. 



3.1 The Model 

Consider a fluid consisting of A'' interacting particles of mass m confined in a two- 
dimensional box of size L x L with periodic (toroidal) boundary conditions. The 
particles interact via a truncated and shifted Lennard- Jones (LJ) pair potential 

[ 0, if r>rc, ^ ' 

where <j)Lj{r) — 4e [(cr/r)^^ — (cr/r)^] is the LJ potential, r is the interparticle 
distance, and Tc is the cut-off which we shall set at Tc = 2.5a. The parameters 
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a and e are, respectively, the characteristic length and energy. For simulations, 
all the quantities were reduced according to e and a, and ks and m are set to 
unity. 

The uniform (in space and time) external driving field E is implemented by 
assuming a preferential hopping in the horizontal direction. This favors particle 
jumps along the field, as it the particles were positively charged; see dynamic 
details in Fig. El As in the lattice counterpart, we consider the large field limit 
E ^ oo. This is the most interesting case because, as the strength of the field is 
increased, one eventually reaches saturation, i.e., particles cannot jump against 
the field. This situation may be formalized by defining the transition probability 
per unit time (rate) as 

uj{r] 77'; E,T)^^[l + t&nh{E ■ S)] ■ min {1, exp{~A<P/T)} . (3) 

Here, any configuration is specified by 77 = {ri, • ■ • , vn}, where is the position 
of the particle i, that can move anywhere in the torus, (l>{r]) = X]i<j 4>{\^i ~ 
stands for the energy of 77, and S — (x\~ Xi) is the displacement corresponding 
to a single MC trial move along the field direction, which generates an increment 
of energy = ^(ji) — ^(r\). The biased hopping which enters in the first term 
of Eq. (O makes the rate asymmetric under rj ^ i]'. Consequently, Eq. (O, 
in the presence of toroidal boundary conditions, violates detailed balance. This 
condition is only recovered in the absence of the driving field. In this limit the 
rate reduces to the Metropolis one, and the system corresponds to the familiar 
truncated and shifted two-dimensional L J fluid |29I3()| . Note that each trial move 
concerning any particle will satisfy that < |r^ — r^j < Smax, where Smax is 
the maximum displacement in the radial direction (fixed at Smax = 0.5 in our 
simulations). 




Fig. 3. Schematic representation of the accessible (shaded) region for a particle trial 
move at equilibrium (left) and out-of-equilibrium (right), assuming the field points 
along the horizontal direction (x). The right hand side shows typical steady state 
configurations above (upper snapshot) and below (lower snapshot) criticality in the 
large field limit. 
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MC simulations using the rate defined in Eq. |(2Jl show highly anisotropic 
states (see Fig. O below a critical point which is determined by the pair of 
values (/Coo, Too)- A linear interface forms between a high density phase and its 
vapor: a single strip with high density extending horizontally along x throughout 
the system separates from a lower density phase (vapor) . The local structure of 
the anisotropic condensate changes from a strictly hexagonal packing of particles 
at low temperature (below T — 0.10), to a polycrystalline-like structure with 
groups of defects and vacancies which show a varied morphology (e.g., at T = 
0.12), to a fluid-like structure (e.g., at T = 0.30,) and, finally, to a disordered 
state as the temperature is increased further. This phenomenology makes our 
model useful for interpreting structural and phase properties of nonequilibrium 
fluids, in contrast with lattice models, which are unsuitable for this purpose. 
Skipping the microscopic structural details, the stationary striped state is similar 
to the one in lattice models, however. 

3.2 Transport Properties 

Regarding the comparison between off-lattice and lattice transport properties, 
the left graph in Fig. ^ shows the net current j as a function of temperature. 
Saturation is only reached at jmax = 4(Smaa:/37r when T — > oo. The current 
approaches its maximal value logarithmically, i.e., slower than the exponential 
behavior predicted by the Arrhenius law. The sudden rising of the current as 
T is increased can be interpreted as a transition from a poor-conductor (low- 
temperature) phase to a rich-conductor (high-temperature) phase, which is rem- 
iniscent of ionic currents ^Sj- This behavior of the current also occurs in the 
DLG. Revealing the persistence of correlations, the current is nonzero for any 
low T, though very small in the solid-like phase. From the temperature depen- 
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Fig. 4. Left graph: Temperature dependence of the net current for the driven LJ fluid. 
Right graph: Transverse-to-the-field current profiles below criticality. The shaded 
(full) line corresponds to the current (velocity) profile of the off-lattice model. For 
comparison we also show the current profile of the DLG with NN interactions (circle- 
dotted line). Since each distribution is symmetric with respect to the system center of 
mass (located here at L/2) we only show their right half parts. 
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dence of j one may estimate the transitions points between the different phases, 
in particular, as the condensed strip changes from sohd to hquid (T « 0.15) and 
finally changes to a fully disordered state (T w 0.31). 

The current is highly sensitive to the anisotropy. The most relevant infor- 
mation is carried by the transverse-to-the-field current profile j± , which shows 
the differences between the two coexisting phases (right graph in Fig.QJ- Above 
criticality, where the system is homogeneous, the current profile is flat on the 
average. Otherwise, the condensed phase shows up a higher current (lower mean 
velocity) than its mirror phase, which shows up a lower current (higher mean 
velocity). Both the transversal current and velocity profiles are shown in Fig. 01 
The current and the density vary in a strongly correlated manner: the high cur- 
rent phase corresponds to the condensed (high density) phase, whereas the low 
current phase corresponds to the vapor (low density) phase. This is expectable 
due to the fact that there are many carriers in the condensed phase which allow 
for higher current than in the vapor phase. However, the mobility of the carriers 
is much larger in the vapor phase. The maximal current occurs in the interface, 
where there is still a considerable amount of carriers but they are less bounded 
than in the particles well inside the bulk and, therefore, the field drives easily 
those particles. This enhanced current effect along the interface is more promi- 
nent in the lattice models (notice the large peak in the current profile in Fig.^. 
Moreover in both lattice cases, DLG and NDLG, there is no difference between 
the current displayed by the coexisting phases because of the particle-hole sym- 
metry. Such a symmetry is derived from the Ising-like Hamiltonian in Eq. |^ 
and it is absent in the off-lattice model. 
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Fig. 5. The temperature-density phase diagram (left graph) was obtained from the 
transversal density profile (right graph) for A'' — 7000, p = 0.35, and different tem- 
peratures. The coexistence curve separates the liquid-vapor region (shaded area) and 
the liquid phase (unshaded area). The diamond represents the critical point, which has 
been estimated using the scaling law and the rectilinear diameter law (as defined in 
the main text). 



8 M. Di'ez-Minguito et al. 

3.3 Critical Properties 

A main issue is tlie (nonequilibrium) liquid-vapor coexistence curve and the as- 
sociated critical behavior. The coexistence curve may be determined from the 
density profile transverse to the field p_L . This is illustrated in Fig. |5| At high 
enough temperature above the critical temperature the local density is roughly 
constant around the mean system density [p = 0.35 in Fig. EJ. As T is lowered, 
the profile accurately describes the striped phase of density p+ which coexists 
with its vapor of density p- [p- < p+). The interface becomes thinner and 
less rough, and p+ increases while p_ decreases, as T is decreased. As an or- 
der parameter for the second order phase transition one may use the difference 
between the coexisting densities p+ — p- . The result of plotting p+ and p- at 
each temperature is shown in Fig. El The same behavior is obtained from the 
transversal current profiles (Fig.^J. It is worth noticing that the estimate of the 
coexisting densities p± is favored by the existence of a linear interface, which is 
simpler here than in equilibrium. This is remarkable because we can therefore 
get closer to the critical point than in equilibrium. 

Lacking a thermodynamic theory for "phase transitions" in non-equilibrium 
liquids, other approaches have to be considered in order to estimate the criti- 
cal parameters. Consider to the rectilinear diameter law [pj^ -\- p-)/2 = p^o + 
''o(21x) ~ T) which is a empirical fit extensively used for fluids in equilibrium. 
This, in principle, has no justification out of equilibrium. However, we found 
that our MC data nicely fit the diameters equation. We use this fact together 
with a universal scaling law p+ — p- = ao{Toc — T)f^ to accurately estimate the 
critical parameters. The simulation data in Fig. [S] thus yields poo = 0.321(5), 
Too = 0.314(1), and (3 = 0.10(8), where the estimated errors in the last digit are 
shown in parentheses. These values are confirmed by the familiar log-log plots. 
Compared to the equilibrium case gni, one has that To /Too « 1.46. This confirms 
the intuitive observation above that the field acts in this system favoring disor- 
der. On the other hand, our estimate for the order-parameter critical exponent 
is fully consistent with both the extremely flat coexistence curve which char- 
acterizes the equilibrium two-dimensional LJ fluids and the equilibrium Ising 
value, /3ising = 1/8 (non-mean-field value). Although the error bar is large, one 
may discard with confidence the DLG value /3dlg ~ 1/3 as well as the mean 
field value. This result is striking because our model seems to have the symme- 
tries and short-range interactions of the DLG. Further understanding for this 
difference will perhaps come from the statistical field theory. 

4 Final Comments 

In summary, we reported MC simulations and field theoretical calculations to 
study the effect of discretization in driven diffusive systems In particular, we 
studied structural, transport, and critical properties on the driven lattice gas and 
related non-equilibrium lattice and off-lattice models. Interestingly, the present 
Lennard-Jones model in which particles are subject to a constant driving field is 
a computationally convenient prototypical model for anisotropic behavior, and 
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reduces to the familiar LJ case for zero field. Otherwise, it exhibits some ar- 
resting behavior, inchiding currents and striped patterns, as many systems in 
nature. We have shown that the additional spatial freedom that our fluid model 
possesses, compared with its lattice counterpart, is likely to matter more than 
suggested by some naive intuition. In fact, it is surprising that its critical behav- 
ior is consistent with the one for the Ising equilibrium model but not with the 
one for the driven lattice gas. The main reason for this disagreement might be 
the particle-hole symmetry violation in the driven Lennard-Jones fluid. How- 
ever, to determine exactly this statement will require further study. It also seems 
to be implied that neither the current nor the inherent anisotropy are the most 
relevant feature (at least regarding criticality) in these driven systems. Indeed, 
the question of what are the most relevant ingredients and symmetries which 
determine unambiguously the universal properties in driven diffusive systems is 
still open. In any case, the above important difference between the lattice and 
the off-lattice cases results most interesting as an unquestionable nonequilibrium 
effect; as it is well known, such microscopic detail is irrelevant to universality 
concerning equilibrium critical phenomena. 

We acknowledge very useful discussions with F. de los Santos and M. A. 
Mufioz, and financial support from MEyC and FEDER (project FIS2005-00791). 
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